Parallel execution of quantum gates in a long linear ion chain 
via Rydberg mode shaping 
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We present a mechanism that permits the parallel execution of multiple quantum gate operations 
within a single long linear ion chain. Our approach is based on large coherent forces that occur when 
ions are electronically excited to long-lived Rydberg states. The presence of Rydberg ions drastically 
affects the vibrational mode structure of the ion crystal giving rise to modes that are spatially 
localized on isolated sub-crystals which can be individually and independently manipulated. We 
theoretically discuss this Rydberg mode shaping in an experimentally realistic setup and illustrate 
its power by analyzing the fidelity of two conditional phase flip gates executed in parallel. Our 
scheme highlights a possible route towards large-scale quantum computing via vibrational mode 
shaping which is controlled on the single ion level. 

PACS numbers: 03.67.Lx, 33.80.Rv, 32.80.Qk 



The ability to execute multiple quantum operations 
in parallel is believed to be a fundamental requirement 
for achieving large-scale quantum computation [TH2]. 
Among the many types of systems being considered for 
the physical implementation of a quantum processor [4], 
trapped ions have attracted much attention for the as- 
tonishingly high degree of experimental control that can 
be gained over their internal and external degrees of free- 
dom [3 . One strategy for achieving parallelism is to build 
many local quantum processors. Current proposals en- 
vision setups where ions are confined in spatially sepa- 
rated wells provided by arrays of microtraps [5] or by 
traps with segmented electric field electrodes [6]. Few 
ions trapped within a given well form one of many lo- 
cal quantum processors that can be operated indepen- 
dently and in parallel [7 . Information can be exchanged 
among different local processors by rearranging the po- 
tential landscape such that previously disconnected ions 
have common vibrational modes. Such rearrangement 
is usually achieved by switching voltages applied to the 
ion trap electrodes [6 . In spite of the availability of mi- 
crostructured arrays it remains a challenge to obtain fast 
switching times and a high spatial resolution of the local 
electric fields that would grant a manipulation of the po- 
tential landscape down to the level of a single ion. This 
places limitations on the transport of ions and the syn- 
chronization of remote quantum processors [6]. 

In this work we introduce a scheme that permits the 
execution of multiple quantum gates in parallel on a long 
linear ion crystal. The method relies on the shaping of 
the vibrational crystal modes through the laser excita- 
tion of selected ions into electronically excited Rydberg 
states. Strong coherent forces acting on these excited 
ions |8j |9] effectively break the long crystal into small 
sub-crystals in the sense that vibrational modes emerge 
which are strongly localized on only a few ions. We illus- 
trate the power of this Rydberg mode shaping by thor- 
oughly analyzing the fidelity of two two-qubit conditional 
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FIG. 1. (Color online) (a) Level structure of Ca + and 
schematics of the envisioned setup. Red and blue symbols 
refer to ions in Rydberg states and ions in electronically low 
lying (ELL) states, respectively. Sub-crystals of ion pairs are 
isolated within a linear crystal formed by 100 ions by the 
excitation of selected ions to the Rydberg nPi/ 2 -state (here 
the 45th/48th and 53rd/56th). Using laser-induced spin- 
dependent forces quantum gates can be executed on the two 
sub-crystals in parallel, (b) Vibrational modes of a crystal 
formed by 100 ions in ELL states. Depicted is the modulus of 
the normal mode matrix B^' 3 ^ where j (m) refer to the mode 
(ion) index (see text for further detail), (c) Vibrational modes 
in the presence of four Rydberg ions. The white dashed lines 
delimit the region corresponding to the ions that are shown 
in panel (a). The Rydberg ions drastically reshape the vibra- 
tional mode structure leading to the emergence of modes that 
are localized on the two sub-crystals (see inset). 



phase flip (CPF) gates that are executed in parallel on 
different sub-crystals belonging to the same ion chain. 
A feature of our scheme is that decoupling between the 
localized modes and the remaining spectator modes per- 
mits us, not only to achieve a high gate fidelity, but also 
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to drastically reduce the complexity of gate optimization 
protocols [T0HT2] . In view of the intrinsic stability of ion 
crystals [13], the ability to address single ions individu- 
ally by lasers [14] and the long lifetime of Rydberg states, 
we believe that this approach indeed highlights a viable 
route towards achieving large-scale quantum computa- 
tion. 

Before providing details let us briefly outline the setup 
we have in mind. We consider a long linear crystal of 
100 40 Ca + ions which is realized within a quartic electric 
potential (see details below). This choice is motivated 
by the proposal discussed in Ref. [13] which envisages 
the implementation of an ion quantum processor where 
a long ion chain is divided into two parts: Quantum com- 
putation is carried out in the central region where ions 
are nearly uniformly spaced. The remaining outer ions 
are continuously Doppler cooled to prevent heating. To 
describe the internal structure of the ions we consider the 
four states depicted in Fig. [IJi. The electronically low- 
lying (ELL) 5, P and D states are employed in numerous 
ion trap experiments for the storage, manipulation and 
read out of quantum information [3]. Furthermore, we 
consider the Rydberg state nPj (with J = 1/2 and the 
principal quantum number n) which is excited from the 
3D 3 /2-state via a single photon transition [T4| [T5]. 

We envisage Rydberg excitations to be carried out in 
the central region of the ion chain and in Fig. |TJl we il- 
lustrate a situation where four Rydberg ions enclose two 
pairs of ions in ELL states. Those ion pairs will form 
the sub-crystals on which we are going to execute quan- 
tum gates in parallel. The underlying physical mecha- 
nism which we aim to exploit for this purpose becomes 
apparent in Figs. []J>,c. Here we show the absolute val- 
ues of the normal mode matrix of the vibrational crystal 
modes which provide a measure on how much each ion 
contributes to a vibrational mode. In Fig. [I}}, that shows 
the case in which all ions are in ELL states, we see that 
in general many ions contribute to each normal mode. 
Compared to this, the presence of Rydberg ions leads to 
a drastic change of the mode structure as can be seen in 
Fig. The reason is rooted in the large polarizability 
V n of Rydberg states [9] which modifies the local trap- 
ping potential, leading essentially to a constriction of the 
ion chain at positions where Rydberg ions are excited. 
The selective Rydberg ion excitation thus creates local- 
ized modes, primarily occupying the two isolated sub- 
crystals composed by ions in ELL states. This Rydberg 
mode shaping permits the parallel execution of quantum 
gates on the two sub-crystals. This is similar in spirit to 
the idea underlying segmented ion traps [I6j [17] . The ad- 
vantage of our approach is that due to the availability of 
single ion laser addressing it fundamentally permits the 
control of the potential landscape on the smallest achiev- 
able length scale - namely on the level of single ions. 
Long linear ion chain - Let us now provide a more de- 
tailed discussion of the practical implementation of the 



above idea. To achieve a long ion crystal we consider 
an ion trap formed by the time-dependent electric po- 
tential $(r,£) = $ r f(r,t) + $ st (r). Here $ r f(r,£) = 
acosQt(x 2 — y 2 ) is the potential of a radio-frequency 
(rf) field with gradient a and frequency Q and 
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x 2 + y 2 , is a quartic static electric potential 



whose parameters f3j (j = 2,4) depend on the specifics 
of the field generating electrodes, i.e. the gradient and 
higher derivatives of the field. Recently similar poten- 
tials have been realized experimentally [I6j [17]. For a 
sufficiently fast rf frequency drive [18] an ion of mass 
M experiences the ponderomotive potential V p (r) = 
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Within this trap an ion chain is 

formed along the z-axis provided that a ^> 1/^2 1, P4I 2 > 0. 
Here l s = [e/ (87reo |/^2 1)] 1//3 is the typical length scale 
associated with V p (r) and e and eo are the elementary 
charge and the vacuum permittivity, respectively. 

Employing a quartic potential has the advantage that 
it allows us to avoid the zigzag instability [T3l [T9J [20] , nor- 
mally encountered in harmonic trapping potentials [21]. 
In addition, the tuning of the parameters fa and fa per- 
mits us to achieve a long ion crystal in which the equilib- 
rium positions of ions in ELL states are approximately 
evenly spaced [13 . The equilibrium positions of the long 
ion crystal is determined by a parameter = 2fal 2 /\fa\. 
In the following we set k± = 1.343 as this choice min- 
imizes fluctuations of the nearest neighbor separation 
within the central region of the ion chain. 

In Ref. [9] we showed that ions excited to the nP^- 
Rydberg state experience not only the ponderomotive po- 
tential but are also subject to an additional radial poten- 
tial that is proportional to the Rydberg polarizability: 
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where T n ~ —0.25 x n 7 (in atomic units). There are 
also small corrections to the trapping potential along the 
z-axis but those are negligible in this linear ion trap. 
The ratio of the radial trap frequencies experienced by 
an ion in the Rydberg/ELL state is approximately given 
by ^Ryd/^ELL = Vl — MVt 2 V n . In practice ratios on the 
order of 2 and larger can be achieved. In the following 
we will show that this is already sufficient for a Rydberg 
ion to effectively introduce a constriction of the linear ion 
chain which strongly affects the vibrational mode struc- 
ture. 

Collective modes and mode shaping - Let us now demon- 
strate the mode shaping considering the transverse 
phonon modes along the x-axis as an example. The 
treatment of the y-phonons is done accordingly. The 
phonon Hamiltonian is [7 H v = X)j=i ^ w j( a ] a j + V^)- 
Here a), (aj) is the creation (annihilation) operator of 
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the j-th phonon, whose frequency uoj is calculated by di- 
agonalizing the the Hessian matrix (J^ m 7-L mn ~Bm X ^ 
(uj/us^B^) with 



0.2 0.4 0.6 



and B^,' x; denoting the eigenvectors. The parameters z m 
are the z-component of the equilibrium position of the m- 

(x) 

th ion and Um is the state-dependent trapping frequency 
of the ra-th ion along the x-axis, i.e. = uj^ydi^ELh) 
if the ion is in the Rydberg (ELL) state. For convenience, 
we have defined a reference frequency uj s = ^J2e\^2\|M 
such that the Hessian is dimensionless. 

We start with a simple situation where in our chain of 
100 ions the 45th and 56th are excited to the Rydberg 
state. The resulting change of the mode structure be- 
comes directly apparent in the modulus of B^,'^ which 
is depicted in Fig. [2ji. Compared to the situation with- 
out mode shaping (Fig.J^)), the striking difference is that 
the 46th to 55th ion constitute a virtually isolated sub- 
crystal hosting a series of spatially localized modes. The 
energies of these local modes are shown in Fig. where 
we also undertake a comparison to the mode energies ob- 
tained by considering exclusively the sub-crystal, i.e., a 
truncated linear crystal composed out of only 10 ions. 
Note, that although we are considering here a case in 
which the sub-crystal ions are symmetrically positioned 
around the center of the ion crystal, our observations re- 
main true also in asymmetric situations. To effectively 
create sub-crystals with localized modes we require that 
(ojRyd/oJs) 2 > max(H mn ) (m ^ n) with max(H mn ) be- 
ing the maximum of the off-diagonal matrix elements of 
the Hessian. This condition means that the energy of 
vibrational modes to which Rydberg ions participate sig- 
nificantly is much larger than the energy of the collective 
modes of ions in ELL states. The mode shaping mecha- 
nism is not limited to the above examples, but also works 
in more intricate situations. An example was given at the 
beginning of this work for the case in which four Rydberg 
ions delimit two sub-crystals. 

Parallel CPF gates - Let us return to this initial exam- 
ple in which we had two sub-crystals composed by the 
ion pairs {46,47} and {54,55}. Each of the sub-crystals 
hosts two localized vibrational modes. The eigenvector 
corresponding to the localized mode with higher energy 
is displayed in the inset of Fig. In the following we 
will show that with these local modes, we can execute 
two two-qubit gates in parallel. Specifically, we discuss a 
a z -type [22] two-qubit CPF gate. Qubits are encoded in 
two ELL states of an ion, denoted by | t) and | \). These 
can be hyperfine states as discussed in Refs. [23, 24], or 
states coupled by optical quadrupole transitions as, e.g., 
in Ref. [25]. The CPF gate is implemented by a laser 
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FIG. 2. (Color online) (a) Modulus of the eigenvector x) 
when ions 45 and 56 are in the Rydberg state. Localized 
modes reside on ions in ELL states forming the sub-crystal 
delimited by the two Rydberg ions. Panel (b) shows the cor- 
responding eigenenergy of the localized modes obtained from 
the full (circle) and truncated (diamond) calculations. The 
largest discrepancy between these calculations is about 0.3%. 
The dots show the (quasi-continuous) energy spectrum of the 
ion chain without mode shaping. In all the calculations we 
use cjell/cJs = 150 and cjR y d/cj s = 198.5. 



induced coupling (see gate lasers in Fig. [TJl) between the 
qubit states and the vibrational crystal modes/phonons. 
This results in a 'spin-dependent' force [23H25] whose ac- 
tion is described by the spin-phonon Hamiltonian [TTJ [26] 

N 

H I= J2 fifim(t)^ ) B«' x >(oJe^ t + h.c). (3) 



Here ft m (t) is the time-dependent Rabi- frequency of the 

(i) 

gate laser that addresses the ra-th ion and rjm = klj 
is the corresponding Lamb-Dicke parameter, with k be- 
ing the modulus of the laser wave vector and lj = 
^h/2Mujj being the oscillator length associated with 
the j-th phonon mode. The CPF gate is conducted by 
switching the gate lasers on for a given time r during 
which the Rabi frequencies fl m (t) are varied. Using the 
Magnus formula [27] , the evolution operator due to Hi is 
then given by 



U (r) = exp 



The first term in the exponential characterizes the resid- 
ual coupling of the ra-th qubit with the phonon modes 
and depends on Q m (r) = ^j[o^m (t)o)- + h.c] where 

(Xrnij) is a parameter that characterizes the coupling 
strength. The second term gives rise to a phonon- 
induced spin-spin coupling between the ra-th and n-th 
qubit thereby effectuating a CPF gate. A perfect CPF 
gate is realized when mn (r) = 7r/8 and a$ (t) = 0. As 
shown in Refs. [10-12 this can be achieved via optimiz- 
ing the time-dependent profile of the Rabi frequencies 
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Q m (t). Such optimization is challenging since in general 
many phonon modes contribute even when only a single 
CPF gate operation is conducted within a long ion chain. 
Rydberg mode shaping has the potential to drastically 
reduce the complexity of such optimization procedure as 
even in long crystals only few vibrational modes actually 
couple to the qubit ions located on a sub-crystal. 

Let us now analyze the performance of two CPF gates 
that are executed in parallel on the two sub-crystals de- 
picted in Fig. [T^l. We use the high energy mode (inset of 
Fig. [TJ:) as 'quantum bus'. To assess the performance of 
the gate operation we consider the following initial con- 
dition: The phonons are in a thermal state, described by 
density matrix p v [18 and the qubits are prepared in a 
product state |tf(0)) = (|^ mi ) <g> |^ ni » <8> (\ip m2 ) <8> |^n 2 )) 
with |VVn) = (| tm> + I im))/V2 and {m^rij} being in- 
dices of ions forming the j-th sub-crystal. Ideally, the 
output state after the parallel execution of the two CPF 
gates is |*(r)> = exp[m/4« i < l + <X a )]|*(0)>. 
However, due to the residual phonon-qubit coupling 
this state will be only reached with a certain proba- 
bility, which we characterize through the fidelity F = 
(*(r)|Ti: v p(r)|*(r)). Here, p(t) = U(t)p(0)W(t) with 
p(0) = p v ® |*(0)) (^(0)| and Tr v denotes trace over the 
vibrational modes. The gate lasers are switched on for a 
time r = 8rb, where Tb = 2tt /uj^ is the oscillation period 
of the bus mode and the Rabi frequency is assumed to 
follow Qrn(t) — sin(z4) (as also discussed in Ref. [24]). 
We optimize the fidelity with respect to the parameter v 
of this simple ansatz. Imposing a two-qubit phase shift 
4>mn(j) = 7r/8 fixes the value of the amplitude (for 
more detail see Ref. [26]). 

Let us first consider a situation in which the two CPF 
gates start simultaneously. We find that the highest 
achievable fidelity within our simple ansatz is F max w 
99.96%. As shown in Fig. [3^i, this maximum occurs at 
vr = 2ir x K (with K an integer). At these points 
the bus modes almost entirely return to their initial 
states [TTJ [28]. The fact that such a high fidelity is 
achievable within this simple ansatz is a direct conse- 
quence of the fact that the Rydberg ions delimiting the 
sub-crystal lead to a dramatic reduction of the number of 
vibrational modes that couple to the qubit ions. Without 
this Rydberg mode shaping the highest fidelity that we 
can achieve is 93%. 

The power of the mode shaping becomes even more ap- 
parent when we relax the requirement of a simultaneous 
start of the two CPF gates. We calculate the fidelity of 
the two CPF gates by introducing a start time delay t& of 
the second CPF gate with respect to the first one. F max 
slightly decreases with growing t& but always remains 
above 98.7% as shown in Fig. [3)3. This demonstrates 
that the two CPF gates can be operated essentially in- 
dependently. In the absence of mode shaping, however, 
F max quickly drops with increasing t& reaching a minimal 
value of « 36.1%. 
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FIG. 3. (Color online) (a) Fidelity of the two CPF gates. The 
solid (dotted) curve is the result calculated using all (only the 
four localized) modes. The dashed curve corresponds to the 
gate fidelity without mode shaping, where the highest energy 
mode (Fig. [Sj)) is used as bus mode. In the calculations, we 
assume that all phonons have the same temperature, with 
the average phonon number of the bus mode being 3.25. (b) 
Maximal fidelity vs. delay time. F max is found by maximizing 
the fidelity over v within the range shown in (a). The solid 
(dashed) curve stands for the calculation with (without) mode 
shaping. 



Finally, let us discuss additional sources that would 
influence the gate fidelity. Firstly, the fidelity will be 
in principle reduced by the radiative decay of the Ry- 
dberg state. This is mitigated by the fact that CPF 
gates are inherently fast [26] and that Rydberg states 
are long lived. For example, choosing the 6OP1/2 state, 
we obtain a Rydberg lifetime « 270 ps and with the trap 
parameters a = 7 x 10 8 V/m 2 , Q = 2ir x 25.2 MHz and 
/?2 = —2.09 x 10 3 V/m 2 we obtain a gate time r ~ 3.7 ps. 
Thus the gate fidelity will be modified by an overall fac- 
tor, about 0.986 iVRyd with A^R y d the number of Rydberg 
ions that are excited during the gate operation. How- 
ever, in principle the gate operations can be accelerated 
by using a more sophisticated optimization protocol [13 , 
and furthermore optimized gate schemes can be imagined 
in which the Rydberg ions do not stay permanently ex- 
cited. Secondly, infidelities are caused by other factors, 
such as the anharmonicity of the ionic motion and cor- 
rections beyond the Lamb-Dicke limit. These have been 
investigated in detail by Lin et al. in Ref. [13^ and their 
contributions have been found to be marginal. 

In conclusion, we showed that the transverse vibra- 
tional modes of a linear ion chain can be shaped by the 
selective excitation of Rydberg ions leading to the emer- 
gence of strongly localized modes. This permits the par- 
allel execution of two- or multi-qubit quantum gates and 
has the potential to drastically reduce the complexity of 
numerical gate optimization schemes for large crystals. 
Rydberg mode shaping can, in principle, be applied also 
to higher dimensional ion crystals and might indicate a 
route towards scalable quantum computation within a 
single large ion crystal. 
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